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Atomic and molecular gas in the merger galaxy NGC 1316 
(Fornax A) and its environment 
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Abstract. We present and interpret observations of atomic and molecular gas toward the southern elliptical galaxy 
NGC 1316 (Fornax A), a strong double-lobe radio source with a disturbed optical morphology that includes 
numerous shells and loops. The 12 CO(1-0), 12 CO(2-l), and Hi observations were made with SEST and the VLA. 
CO emission corresponding to a total molecular hydrogen mass of ~5xl0 8 Mq was detected toward the central 
position as well as northwest and southeast of the nucleus in the regions of the dust patches. The origin of that 
gas is likely external and due to accretion of one or several small gas-rich galaxies. Hi was not detected in the 
central region of NGC 1316, but ~2xl0 7 M© of atomic gas was found toward the giant Hn region discovered by 
Schweizer (1980) located 6'7 (or 36.2 kpc) from the nucleus. Hi was also found at three other locations in the 
outer part of NGC 1316. The Hi distributions and kinematics of the two nearby spiral companions of NGC 1316, 
NGC 1317 (a barred galaxy to the north) and NGC 1310 (to the west) could be studied. Both galaxies have 
unusually small Hi disks that may have been affected by ram-pressure stripping. 

Key words, galaxies: ellipticals; individual (NGC 1316, NGC 1317, NGC 1310); ISM; kinematics; interactions- 
galaxies: clusters: Fornax cluster- radio lines: galaxies 



1. Introduction 

Over the past fifteen years, evidence for the existence of 
a complex, multi-phase interstellar medium in elliptical 
galaxies has accumulated. Most of the gas in ellipticals is 
hot (10 7 K) and radiates at X-ray wavelengths, but there 
is also a cooler component of Hi , CO, ionized gas and dust 



(e.g. Bregman et al. 1992). The origin and fate of the gas 
are still debated, but it seems that the hot gas is associated 
with the spheroid, whereas the atomic and molecular gas 
as traced by the CO characterize a prominent disk (e.g. 
van Gorkom & Schiminovich 1997). The detection rates 



of Hi or CO are significantly higher among ellipticals with 
peculiar features (such as tidal tails, loops, shells, dust 
lanes) than in the whole class of ellipticals, which suggests 
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that the gas is of external origin and due to accretion of 
smaller companions and/or to major mergers. 

NGC 1316 (Fornax A, Arp 154, PKS 0320-37) on the 
outskirts of the Fornax cluster presents several character- 
istics of a merger remnant. The disturbed outer morphol- 
ogy with numerous shells and loops was first noted by 
Murdin et al. ( |1979| ). The galaxy was extensively stud- 
ied in the optical by Schweizer ( |1980| , hereafter S80) and 
Schweizer ( 198l[ ) who classified it as a D-type galaxy with 
an elliptical-like spheroid surrounded by a large envelope. 
The brightness profile approximately obeys the r 1 ' 4 law. 
The galaxy has an inner dust lane of length ~2' oriented 
close to the minor axis. S80 also reported the discovery 
of a giant Hn region (hereafter denoted SH2) 6.'7 from 
the nucleus, and the existence of an inner disk of ionized 
gas probably associated with the dust lane and extend- 
ing out to 54". The disk is in rapid rotation (wsini ~ 350 
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Fig. 1. Radio continuum contour map at A20 cm of Fomalont et al. (1989ft superimposed on an optical image 
taken from the Digitized Sky Survey. The resolution of the radio continuum map is 14". The two neighbor galaxies, 
NGC 1317 and NGC 1310, and the location of the extended emission line region EELR and the giant Hn region SH2 
in the southern tidal loop called Li by Schweizer ( 198CJ ) are indicated. NGC 1310 is located in the foreground of the 
western radio lobe. The black rectangle shows the size of the image displayed in Fig. 2. 



km s _1 ). To explain those features and the lack of a nearby 
perturber, S80 suggested that NGC 1316 accreted one or 
several gas- rich companion gala xies a bout 4x 10 8 - 2 x 10° 
years ago. Mackie & Fabbiano ( 1998 ) reported the discov- 
ery of an extended emission-line region (hereafter denoted 
EELR) located in a tidal tail 5f3 from the nucleus and sug- 
gested on the basis of its position and extreme size that it 
is related to the merger process. 

NGC 1316 is the third nearest strong radio galaxy 
in the sky after NGC 5128 (Centaurus A) and M 87 
(Virgo A) (see Fig. 1). The radio emission arises in two 
giant lobes separated by 33' and located mainly outside 
the o ptical galax y, at a position angle of ~ 110° (Wade 



196l| ; Ekers et al. |l983j , Fomalont et al. |l989|) . In the cen- 



tral arcminute a bent radio jet abruptly terminates at 
the location of a dust la ne, wh ich may have deflected it 
(Geldzahler & Fomalont 1984 ). The nucleus is a br ight 
UV source, unresolved by HST (Fabbiano et al. 1994 ). 



A study of the stellar kinematics of NGC 1316 showed 
that the galaxy is rotating rapidly around the minor axis 
at a position angle of ~ 140° (Bosma et al. 1985 ). More 
recently Arnaboldi et al. ( 1998| ) measured radial velocities 
of 43 planetary nebulae in the outer parts of NGC 1316 
to determine the distribution of angular momentum. They 
derived a constant velocity v= 110 km s _1 along the major 
axis of the galaxy in the radial range i?=50-200". They 
also probed the kinematics of the inner region (R < 80") 
using absorption line spectra. They confirmed the align- 
ment between the photometric and the kinematic minor 
axes but also found evidence for a disturbed kinematical 
pattern beyond R ~ 30". 

X-ray emission from NGC 1316 has been observed with 
several instruments. The Lx/Lb ratio of NGC 1316 is low 
compared to that of E and SO ga laxies {Eins tein o bser- 



vations by Fabbiano et al. 1992| ; Kim et al. 1992J) , but 



in the same range as that of dynamically young ellipti- 
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cals which may have experienced mergers (Fabbiano & 
Schweizer 1995; see also Sansom et al. 2000). 

In order to place NGC 1316 into the merger sequence 



Table 1. VLA observing parameters 



as discussed by Hibbard & van Gorkom (1996) it is nec- 



essary to know the distribution of the cool gas, both 
the neutral atomic hydrogen and the m olecul ar hydro- 
gen as traced by the CO. Sage & Galleta ( 1993 ) detected 
12 CO(1-0) line emission at two positions along the inner 
dust lane and with kinematics similar to that of the ionized 
gas. In this paper, we present more extensive and more 
sensitive maps of the 12 CO(1-0) and the first 12 CO(2-l) 
observations in the central part of NGC 1316 as well as a 
search for CO emission toward the outer Hil complexes. 
We also present results from a VLA search for Hi gas in 
NGC 1316 and its surroundings. 

The paper is organized as follows: We first describe 
the Hi and CO observations, then discuss the Hi results 
for the companion galaxies NGC 1317 and NGC 1310. We 
then focus on the Hi and CO results for NGC 1316 and 
the implications for the merger picture. 



2. Observations and data reduction 

2.1. Hi observations 

The Hi observations consist of data collected with the Very 
Large Array (VLA)F] in its spectral-line mode. The galaxy 
was observed in the hybrid DnC configuration: the 1.3 km 
D-array configuration with the antennas of the north arm 
in the 3 km C-array configuration to improve the North- 
South angular resolution, thus compensating for the low 
declination of NGC 1316. The observations consisted of 
two runs of 5 hours each. 

A 6.25 MHz bandwidth was used, corresponding to 
a total velocity range of 1300 km s — 1 . On-line Hanning 
smoothing has been applied to the data after which ev- 
ery other channel was discarded, leaving a set of 31 inde- 
pendent channels centered at 1774 km s _1 and spaced by 
41.7 km s _1 . The observing parameters are summarized 
in Table 1. 

The data were inspected, calibrated, and maps were 
produced using NRAO's Astronomical Image Processing 
System (AIPS). The continuum was estimated from chan- 
nels free of line emission and was subtracted in the 
UV plane. Line channel maps were made using uniform 
weighting with a ROBUST factor of 1 which, after some 
experimentation, was found to give the best combination 
of sensitivity and spatial resolution. The channel maps 
were CLEANed and a correction for primary beam at- 
tenuation was applied to the entire data cube. Maps of 
the atomic hydrogen column density and velocity field 
were produced. The maps were not smoothed. Only pixels 
above a flux cut-off of 1.5 mJy/beam (2.5 a) were included 
in the summation of the channel maps. 



Phase center (B1950) 
Velocity center 
Primary beam (FWHM) 
Flux calibrator 
Phase calibrator 
Bandpass calibrator 
Array 
Date 
Synthesized beam 

- FWHM: major x minor axis 

- Position angle 
Bandwidth 
Number of channels 
Channel separation 
Time on source 
Noise level (lcr)* 

- Flux density 

- Column density 

- Hi mass 



03 h 20 m 47:2; -37°23'08'.'2 

1774 km s _1 (heliocentric) 

30' 

3C 48 

0332-403 

3C 48 

DnC 

19, 21 Oct 1985 

52" x 41" 

+2° 

6.25 MHz 

31 

41.7 km s" 1 

6 hours 51 min 

0.6 mjy beam _1 ch _1 
1.3xl0 19 cm -2 beam^di -1 
2xl0 6 M beam^di -1 



* Noise level in the continuum-subtracted maps before primary 
beam correction 



The data were taken in 1985 by van Gorkom, Ekers, 
Schweizer and Wrobel. Due to limitations in on-line and 
off-line computing power only 25 antennas and one po- 
larization (RR) could be used. The data were fully re- 
duced by van Gorkom, but never published. Although 
NGC 1310, NGC 1317 and the Hil region were unam- 
biguously detected, the quality of the data reduction was 
poor. This paper shows that reprocessing of ancient (but 
venerable) data is very worthwhile, especially when a large 
spectral dynamic range is needed, because there have been 
major improvements in software such as the UVLIN pro- 
gram within AIPS and far better weighting schemes in the 
imaging tasks. 

2.2. CO observations 

The observations were done with the 15m Swedish-ESO 
Submillimeter Telescope (SESTf] on La Silla, Chile, in 
two sessions on 1999 December 7-13 and in 2001 January. 
We used the IRAM dual receiver to observe simultane- 
ously at the frequencies of the 12 CO(1-0) and 12 CO(2-l) 
lines. At 115 GHz the half-power beamwidth of the tele- 
scope is 43", which corresponds to a linear scale of about 
3.9 kpc at the adopted distance of 18.6 Mpc. The CO(2-l) 
map has a resolution of 22" (or ~2 kpc). This is the high- 
est resolution map of molecular gas in NGC 1316 available 
to date. A balanced on-off beam-switching mode was used. 
The pointing was regularly checked on nearby SiO masers. 
The temperature scale used throughout this article is the 
main-beam temperature T m b = T\jr\ m \,. The main-beam 



1 The VLA is a facility of the U.S. National Radio Astronomy 
Observatory, which is operated by Associated Universities Inc. 
under contract with the U.S. National Science Foundation 



The SEST is operated jointly by the European Southern 
Observatory (ESO) and the Swedish National Facility for 
Radio Astronomy, Chalmers Centre for Astrophysics and Space 
Science 
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Table 2. Basic parameters of the galaxies and outer Hn regions 





Note 


NGC 1316 


NGC 1317 


NGC 1310 


SH2 




EELR 


RA(1950) 


(1) 


03 h 20 m 47.2 s 


03 h 20 m 49.8 s 


03 h 19 m 08.8 s 


03 h 20 m 43.5 s 


03 h 20 m 34.8 s 


DEC(1950) 


(1) 


-37° 23' 08" 


-37°16'50.9" 


-37°16'42.3" 


-37°29'28" 


-37° 27' 39" 


Wheiio [km s" 1 ] 


(2) 


1783 


1948 


1739 


1690 




D [Mpc] 


(3) 


18.6 


18.6 


18.6 


18.6 


18.6 


Type 


(4) 


.LXSOP. 


.SXR1.. 


.SAS5*. 






Size 


(5) 


12' x 8.'5 


2^75 x 2.'4 


2^0 x 1.'55 


1'.'5 x 3'.'5 


81" x 27" 


P.A. [°] 


(5) 


50 


78 


95 




~45 


B° T [mag] 


(6) 


9.4 


11.81 


12.56 






Lb [L q ] 


(6) 


9.35xl0 10 


1.02xl0 10 


5.1xl0 9 






LpiR [Lq] 


(7) 


2.05xl0 9 


2.45 xlO 9 


7.66x10 s 






M(Hi) [M ] 


(8) 


<10 8 


2.7x10* 


4.8x10 s 


1.9xl0 7 


<3xl0 6 


M(Ha) [M ] 


(9) 


~5xl0 8 


3.1x10 s 


<0.2xl0 8 


<7xl0 6 


<7xl0 6 


whijAuhi [km s" 


'] (8) 




1960, 132 


1827, 144 


1690, < 42 


1750, < 42 


Fm [Jykms- 1 ] 


(8) 


<0.1 


3.3±0.5 


5.9 ± 1.3 


0.2 


0.1 


i>co, A«co [km s 


- 1 ] (9) 




1948, 84 








Notes: (1) NGC 


1316: Sc 
tie offsets 


hweizer ( 


1981 


); NGC 1317 
e central posi 


da Costa et 
Lion given by 


al. (1998); NGC 1310: Lo^ 


calculated from t 


relative to th 


Mackie & Fabbiano ( 


1998). 



(2) NGC 1316: Arnaboldi et al. (1998); NGC 1317: Horellou et al. (1995); NGC 1310: RC3; SH2: this work. 



Loveday (1996); SH2 and EELR: 



1999). At that 



(3) Distance estimate for NGC 1365 based on HST measurements of cepheid variable stars (Madore et al. 
distance, 1' corresponds to 5.4 kpc. That distance was adopted for all the systems discussed here. 

(4) RC3. 

(5) For the three gal axies: optical diameters at the 25 mag arcsec -2 isophote taken from RC3; SH2: Schweizer (1980); EELR: 
Mackie & Fabbiano J1998J ). 

(6) RC3. 

(7) Derived from the fluxes at 60 and 100 /xm listed in the IRAS catalogue using the relation: Lfir — 3.94 x 10 5 D 2 [2.58/60 +/100] 
where D is the distance in Mpc, /60 and /100 are the fluxes in Jy, and the far-infrared luminosity Lfir is in solar luminosities. 

(8) This work. 

(9) NGC 1316, SH2 and EELR: this work. H 2 mass of NGC 1317 and NGC 1310: Horellou et al. fll995| ) after correction for the 
distance adopted here. 



efficiency of SEST is r) mh = 0.7 at 115 GHz and 0.5 at 230 
GHz. The data were analyzed using the software CLASS. 
The spectra were smoothed to a final velocity resolution of 
~14.5 km s _1 . Only first-order baselines were subtracted. 

2.3. Determination of gas masses 

The integrated Hi masses are calculated using: 



where 9 is the half-power beam width of the telescope 
and I m b=J T m bdv in K km s _1 . To be consistent with 
previous papers, we have used a conversion factor X = 
N(H 2 )/I(CO) of 2 .3 x 10 20 mol cur 2 (K km s" 1 )- 1 
(Strong et al. 198S). However, use of the Galactic CO- 



H2 conversion factor may lead to an underestimate of the 
molecular gas mass for elliptical gala xies w ith a significant 
amount of cold dust (Wiklind et al. 1995 ). 



M ffl [M G 



2.356 x lO 5 ^ 



Smdv 



where Z?m pc is the distance in Mpc and J Smdv is the 
integrated Hi flux in Jy km s . 

Hi column densities can be derived from the observed 
integrated intensities by using the relation: 



N m = 



110.4 x 10 3 



ah 



Smdv 



where A^hi is the column density in 10 cm per beam 
per channel, and a and b are the diameters of the synthe- 
sized beam at full width half maximum in arcseconds. 

From the 12 CO (1-0) line intensities we can estimate 
the mass of molecular gas (H2) within the 43" SEST beam: 

M(H 2 )[M Q ] = 1.25 x 10 5 (0/43") 2 I mb {CO) D 2 Mpc 



3. Atomic hydrogen 

Figure 2 shows the Hi field of interest around NGC 1316 
superimposed on an optical image from the Digitized 
Sky Survey. Hi emission was found in the two companion 
galaxies NGC 1317 (to the north) and NGC 1310 (to the 
west). No Hi was detected in the main body of NGC 1316 
to a limiting mass of ~ 10 8 M Q but some was found at 
four locations in the outer part of the galaxy. One of these 
is associated with the outer Hll region SH2. The other is 
coincident with the extended emission line region EELR. 
The two other Hi concentrations to the west and north- 
west of NGC 1316 are located at the edge of an optical 
plateau between the area of X-r ay em ission as observed 
by ROSAT (Mackie & Fabbiano |l998[ ) and the extended 
western radio lobe (see Figure 1). The mass of each of 
these Hi concentrations is ~ 2 x 10' Mr. We will defer 
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Fig. 2. Integrated Hi map overlaid on an optical image from the Digitized Sky Survey. The contour levels are n x 
5.2 x 10 19 cm -2 . The first contour is at the 4ct level. The two giant Hn regions to the south, SH2 and EELR, are 
detected. No primary beam correction has been applied to this map. 



discussion of these Hi features and the lack of Hi emission 
in the main spheroid of NGC 1316 until later and first 
discuss the Hi properties of NGC 1317 and NGC 1310. 

3.1. The neighbor galaxies NGC 1317 and NGC 1310 

NGC 1317 is a small barred galaxy north of NGC 1316. 
The two galaxies have a projected separation on the sky 
of 6'3 (34 kpc) and a radial velocity difference of +165 
km s _1 . An inner bar perpendicular to the primary bar 



is visible in the opti cal images of Schweizer ( 198CJ ) and 
Wozniak et al. ( 1995 ). A ring of star formation is seen in 
the optical images and in the ultraviolet image obtained 
with the Ultraviolet Imaging Telescope. The primary bar 
is "fat" (~60"x80"). The secondary bar has a position 
angle of ~60° and a total length of ~14". S80 assigns 
NGC 1317 a late SB(rs)a or early SB(rs)b type in the 
Hubble classification. 

Figure 3 presents in form of isocontours the atomic 
hydrogen distribution and the velocity field in NGC 1317 
derived from the observations. Also shown is the global 
Hi spectrum. The channel maps are shown in Figure 4. 
The Hi distribution is more extended than the beam size 
and is elongated in the direction of the minor axis of the 
galaxy. The gas is clearly rotating, with the south-west 
side receding. The velocity and the linewidth of the global 
Hi spectrum agree with t hose of the 12 CO(1-0) line de- 
tected by Horellou et al. (|l995|) (see Table 2). NGC 1317 
has comparable amounts of molecular and atomic gas 
(~3xl0 8 M ). 



The Hi disk of NGC 1317 is small compared to the op- 
tical extent. The Hiisophote corresponding to a column 
density of 10 20 cm -2 crosses the major axis of the galaxy 
at a radius of about 55" whereas the optical isophote at 
25 mag arcsec -2 is found at 82'.'5. The ratio of the ex- 
tent of the Hi disk to that of the optical image is therefore 
0.7, which places NGC 1317, if classified as an Sab galaxy, 
in the group of galaxies with an Hi distribution strongly 
affected by the environment (group ill of Cayatte et al. 



1994 for Virgo galaxies). Cayatte et al. (1994) showed that 



those Virgo galaxies with a truncated Hi disk have a nor- 
mal central Hi surface density average as compared to field 
galaxies and concluded that the galaxies (all located close 
to the cluster core) have been affected by ram-pressure 
sweeping. NGC 1317 has a central Hi surface density aver- 
aged over half the optical radius of ~3xl0 20 cm -2 , which 
is lower than that of Sab galaxies in the comparison sample 
of Cayatte et al. (5.0±1.2xl0 20 cm -2 ), but in the range 
of values found for galaxies of earlier types (SO, SO/a, Sa 
have 3.4 ±0.9 x 10 20 cm -2 ). Given the uncertainty in the 
morphological type and the fact that we are discussing an 
individual galaxy, it is difficult to draw any firm conclusion 
about the central Hi density of NGC 1317. Nevertheless, 
it is clear that the Hi disk is small and that NGC 1317 has 
a low Hi mass for its optical size and morphological type. 
We calculate an Hi deficiency of .64 in the log, or a factor 
of 4.4 using the average value for isolated Sab galaxies ob- 



served by Haynes & Giovanelli (1984) and the definition: 
HIdef = \og{M(HI)/D 2 ) 



NGC1317 



\og(M(HI)/D 2 )cs 
where cs stands for control sample. Although NGC 1316 
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Fig. 3. a. Hi integrated intensity map of NGC 1317 overlaid on a grey scale optical image from the Digitized Sky 
Survey. A correction for the primary beam has been applied. The levels are (4, 8, 16, 31, 47)xl0 19 cm -2 . The first 
contour corresponds to the 2a level, b. Hi velocity field overlaid on a grey scale image of the Hi integrated intensity. 
Contours of constant heliocentric radial velocity begin at v — 1910 km s _1 in the northeast and increase in steps of 20 
km s _1 toward the southwest. The cross marks the position of the optical center of NGC 1317. c. Integrated Hi line 
profile derived from the observations. 
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Fig. 4. The Hi line emission of NGC 1317 displayed as channel maps running from 2066 km s _1 (top left panel) to 
1857 km s _1 (bottom right panel). North is up and east is left. The levels are 1.5, 3, 4, 5, 6, 7, 8, 9 mJy/beam. The 
first contour corresponds to the 2a level. The cross marks the position of the optical center of NGC 1317. 



is not located in the center of the Fornax cluster, it does 
dominate its corner of the cluster and is a source of X-ray 
emission. It is possible that NGC 1317 has passed through 
the X-ray halo of NGC 1316 and been stripped of some of 
its atomic gas. 

NGC 1310 is an almost face-on spiral fortuitously posi- 
tioned between us and the western radio lobe of Fornax A. 
Ionized gas from NGC 1310 produces a Faraday screen 



that depolarizes the background radiation (Fomalont et al. 
1989; Schulman & Fomalont 1992). Although the galaxy is 



located 20'5 away from the center of Fornax A on the plane 
of the sky and in the outer region of the primary beam 
of the VLA, Hi emission is well detected. Figures 5 and 6 
present the results of the Hi observations in the same form 
as those for NGC 1317. The approaching gas is located on 
the eastern side of the galaxy. If the spiral arms are trail- 
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Fig. 5. a. Hi integrated intensity map of NGC 1310 overlaid on a grey scale optical image from the Digitized Sky 
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Integrated Hi line profile derived from the observations. 
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Fig. 6. The Hi line emission of NGC 1310 displayed as channel maps running from 1857 km s _1 (left panel) to 1732 
km s _1 (right panel). North is up and east is left. The levels are 4, 8 and 12 mJy/beam. The first contour corresponds 
to the 2ct level. The cross marks the position of the optical center of NGC 1310. 



ing, then the northern part of NGC 1310 is the near side. 
The Hi distribution is slightly asymmetric and the emis- 
sion is stronger on the eastern side. We estimate a total 
mass of atomic hydrogen of 4.8 xlO 8 M Q for NGC 1310, 
which is about 20 times higher than the upper limit on 
the H2 mass set by Horellou et al. ( 1995 ). 

Comparing the extent of the Hi disk to that of the opti- 
cal galaxy, one finds that the isocontour of column density 
10 20 cm~ 2 roughly coincides with the blue isophote at the 
25th mag arcsec -2 . This is unusual compared to Sc galax- 
ies in the field, which normally have an Hi diameter larger 
than their D25 diameter by a factor of 1.88±0.08 (Cayatte 
et al. |l994j) . However, NGC 1310 is not significantly Hl- 



deficient comp ared to the isolated Sc galaxies of Haynes 
& Giovanelli ( |l984D (HIdef=0.17 in the log). The small 
size of the disk (~2' across) together with the location of 



NGC 1310 in the outer part of the primary beam of the 
VLA make it difficult to say more about its Hi properties. 

3.2. NGC 1316 (Fornax A) 



• The spheroid. Huchtmeier & Richter ( |1989| ) give an 
upper limit on the Hi integrated line intensity of 4.5 
Jy km s _1 , which corresponds to 3.7xl0 8 M Q at our 
adopted distance. We didn't detect any Hi in the central 
region of the galaxy down to a 5 a limit of 10 7 M Q per 
beam and per channel. 

If, however, the Hi emission is as broad in velocity as 
the CO detected in the central region (see next section) 
then the limit is 8.5 x 10 7 M Q . This could be even higher 
if the Hi is extended over several beams. Very extended 
Hi , moreover, will be missed in the present observations 
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Fig. 7. a. Contour map of the integrated Hi line emission in the Hn region SH2. The levels are 2,3,4 and 5 mJy/beam. 
The first contour corresponds to the 2a level. The circle shows the size of the beam. The cross marks the optical position. 
b. Integrated Hi spectrum of the Hn region SH2. 



because the quality of the data at the shortest spacings is 
limited by the presence of 34 Jy of continuum emission. A 
safe limit to the Hi emission in the main body of NGC 1316 
is 10 8 M . 

We did not detect Hi absorption towards the nuclear 
continuum source which has a flux density of 220 mJy 
at 1412 MHz. The 5 a column density limit is 1 x 10 18 
T s cm -2 , where T s is the spin temperature (in K) of the 
absorbing atomic hydrogen. 

The four Hi detections in the outer regions of 
NGC 1316 are all far outside the main body of NGC 1316 
and lie at or close to the edge of the faint optical shells 
and X-ray emission of NGC 1316. Each clump contains 1 
to 2 xlO 7 Mq of Hi . The southern clumps are associated 
with optical emission line gas: the Hn region SH2 and the 
EELR. 

• SH2. This Hn region has a size of 1"5 x 3"5 (135 pc 
x316 pc); it is located 6f7, or 36.2 kpc south of the nucleus 
and has a radial velocity of -105±8 km s _1 relative to it 
(S80). S80 notes that a "blue fuzz" , probably an associated 
complex of OB stars, is visible near the Hn region. 

The Hi emission in SH2 appears in one channel cen- 
tered at u=1690 km s _1 (our channel width is 41.7 km s _1 ) 
(see Figure 7). This is in agreement with the velocity of 
the ionized gas quoted by Schwcizcr. The Hi mass is about 
2x 10 7 M©. The distribution is unresolved at our resolution 
of 52" x 41" (4.7 x 3.7 kpc). 

• The EELR. This extended emission line region is lo- 
cated 5'3, or 28.6 kpc from the nucleus and has a size of 
81" x 27", or 7.3 kpc x 2.4 kpc (for our distance of 18.6 
Mpc) . On the plane of the sky, it appears on a region of 
high surface brightness in the loop called Li by S80 and 
near a major dust lane. 



The Hi associated with the EELR is - 1 x 10 7 M Q and 
also unresolved by our beam. 

The two northern clumps are located at the outer edge 
of a plateau in the optical surface brightness distribu- 
tion (Fig. 2, see also Fig. 7 of S80). The velocities of the 
Hi clumps are (from north to south through west) 1857, 
2000, 1750 (EELR) and 1690 (SH2), respectively, follow- 
ing the general rotation pattern on the central CO (next 
section) and suggesting that the clumps may be the left- 
overs of a more widespread Hi structure, possibly a tidal 
shred. The location and velocity structure of the Hi are 
reminiscent of other shell galaxies such as Centaurus A 
(Schiminoyich et al. |1994| ) and NGC 2865 (Schiminovich 
et al. |1995| ). 



4. Molecular gas 

4.1. The central region of NGC 1316 (Fornax A) 

Results from the CO observations are presented in 
Figures 8 to 10 and the line parameters are given in 
Table 3. The observations were done by two observers 
in two separate sessions, and the observing strategy was 
slightly different. During the first observing session (see 
top panel of Fig. 8), we mapped the CO emission along 
the axis with a position angle PA=142°, given by S80 as 
the direction of the maximum extent of the ionized gas 
and probably the rotation axis of the spheroid (perpen- 
dicular to the optical major axis). During the second ses- 
sion (Fig. 8, bottom panel), the grid was in RA/DEC. 
Emission is detected north-west and south-east of the nu- 
cleus as well as toward the center. The intensities in the 
NW region are about twice as strong as in the SE, both 
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Fig. 8. Top: CO(l-O) and CO(2-l) spectra obtained in December 1999. The offsets arc in arcseconds. The x scale 
ranges from 1000 to 2300 km s _1 . The y scale (main-beam temperature, T m b) ranges from -0.01 to 0.03 K. The velocity 
resolution of the CO(1-0) spectra is 14.5 km s -1 , that of the CO(2-l) spectra is 14.6 km s _1 . Bottom: Results of the 
January 2001 observing sesion. The x scale ranges from 1100 to 2300 km s _1 . The y scale and the velocity resolution 
are the same as for the spectra above. 

in CO(1-0) and CO(2-l). The CO distribution roughly CO(l-O) spectrum is about 5 times lower t han th at in the 
follows that of the dust along the minor axis and to the spectrum obtained by Wiklind & Henkel ( |1989[ ) and the 
NE. bandwidth is larger, which makes it easier to see the line 

structure. The integrated line intensity we infer is about 



• Emission was detected toward the central position. 
The averaged central spectra obtained during both observ- 
ing sessions are shown in Figure 9. The noise level in our 



three times lower than that reported in that work, but it 
may agree within the not quoted error bars of Wiklind 
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Table 3. Line parameters 
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(1), (2): Offset positions after rotation of the coordinate system by PA=142°. 

(3), (4): Offset positions in RA and DEC. 

(5) to (7): For all positions but the center we give the parameters of the gaussian lines fitted to the CO (1-0) spectra: central 
velocity, linewidth, integrated main-beam intensity. For the center position where the lines are clearly not gaussian, we give 
the velocity of the centroid, the linewidth and the main-beam intensity integrated over the velocity range of the line (see 
Figure 9). For non-detections, the limit on the integrated intensity was taken as 3xa x 150 km s _1 . 

(8): Noise per 14.5 km s _1 channel. 

(9) to (12): Same as columns (5) to (8) but for the CO(2-l) spectra. 



& Henkel. Sage & Galletta ( |i"993| ) had detected CO(1-0) 
at two positions 30" to the SE and 50" to the NW of the 
nucleus but not toward the center. The fluxes and the lim- 
its quoted agree with ours within the uncertainties. The 
CO(1-0) line seems to be very broad (~500 km s _1 at the 
base) and to consist of three components. It is likely that 



the molecular gas distribution is discrete and that the low- 
velocity component arises from emission from the SE that 
enters the 43" beam, whereas the hint of emission at high 
velocity (~1900 km s -1 ) comes from the NW concentra- 
tion. All those features are, if present at all, much weaker 
in the CO(2-l) spectrum taken with a narrower beam. 
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We calculate a mass of the H2 gas of ~10 8 M Q within the 
inner 43" (3.9 kpc). This estimate of mass is based upon 
a standard Galactic conversion factor. 

• The NW concentration is centered at Aa=-14", 
A<5=+37" relative to the central position, that is, at a 
projected distance of 3.5 kpc from the nucleus. The emis- 
sion is extended on the scale of the 22" CO (2-1) beam, 
and the measured apparent diameter at half maximum in- 
tensity is 35", which correspond to a true extent of somce 
= 27", or a linear size of ~2.5 kpc. This is consistent with 
the CO (1-0) measurements, which yield an apparent di- 
ameter of the half peak intensity of 50" (and a true source 
size of 24"). Velocities around 1850 km s _1 are found in 
the NW region, that is, ~70 km s _1 higher than the sys- 
temic velocity of NGC 1316. A slight east-west velocity 
gradient is observed, with higher velocities to the west. 

From the apparent extent of the northwest CO concen- 
tration measured above and the CO(l-O) peak intensity, 
we derive a mass of molecular hydrogen within that region 
of 2.2xl0 8 M Q . 

• In the southeast, the strongest emission is detected 
22" (^2 kpc) away from the nucleus and is half as strong 
at 44". It corresponds to a H 2 mass of 0.86x10 s M Q . The 
velocities measured there are of the order of 1550 km s _1 , 
that is, ~230 km s _1 lower than the systemic velocity of 
NGC 1316. 

• The CO velocity fields are displayed in Fig. 10c and 
lOd. They are fairly regular. The systemic velocity of 
NGC 1316, however, is likely to be higher than suggested 
from those interpolated contour maps. As discussed above 
for the central position, the low velocity wing around 1500 
km s _1 could be due to emission from the NW concentra- 
tion that enters the beam. The centroid velocity of 1559 
km s _1 given in Table 2 for the central position would 
therefore be too low and affect the interpolation in the 
contour map. 



4.2. The outer emission-line regions 

We have searched for CO emission toward the giant Hn 
region SH2 and the extended emission line region but 
have not detected any, down to the noise levels quoted 
in Table 3. For SH2 as for EELR, this corresponds to a 5a 
limit on the H2 mass of ~10 7 M Q per 43" beam per 14.5 
km s _1 channel. 



5. Discussion 

While CO is clearly detected in the central region of 
NGC 1316, atomic hydrogen is not. The distribution of 
the molecular gas follows that of the dust, and its kine- 
matics provides new insight into the dynamical history of 
the system. Hi is detected in the outer part where no CO 
is found down to a similar mass limit. 



5.1. Comparison with other components of the ISM 

• The central part of NGC 1316 
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Fig. 9. CO(l-O) detection toward the center of NGC 1316. 
The line is broad and seems to consist of three compo- 
nents. The y axis shows the main-beam temperature in 
Kelvin. The width of the window set to subtract a base- 
line is shown on the x axis: 1444 to 2032 km s _1 for the 
CO(1-0) spectrum as indicated by the box, and 1460 to 
1917 km s _1 for the CO(2-l) spectrum as indicated by 
the ticks. The CO(2-l) line seems to be narrower than 
the CO(1-0) one, which is to be expected if the molecular 
gas is located in the dusty regions offset from the cen- 
ter and located outside the narrower CO(2-l) beam. The 
spectra have been produced by averaging the data taken 
during both observing sessions. 



The dust patches are very visible in the B — I HST 
image of NGC 1316 displayed in Figure lie. They ap- 
pear on the NW and the SE side of the nucleus, roughly 
aligned along the minor axis of the galaxy. The SE lane 
can be traced to wi thin 1 arcsec from the center, which 
led Schweizer ( 1980J ) to suggest that it lies on the near 
side of the galaxy. The molecular gas distribution at 22" 
resolution follows roughly the distribution of the dust. 

Ionized gas was found by S80 to be associated with 
the dust and distributed in a rotating disk with extremal 
velocities ±350 km s" 1 at 84" NW and 35" SE. The SE 
side of the disk is approaching. The extreme range of CO 
velocities that we measure on either side of the nucleus, 
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Fig. 10. a. and c. Contour map of the CO(l-O) distribution and corresponding velocity field. The levels of the total 
intensity map go from 0.5 to 6 K km s _1 in steps of 1 K km s _1 . The levels of the velocity map go from 1500 to 2000 
km s _1 in steps of 50 km s _1 . The circle shows the size of the beam. The black dots indicate the observed positions. 
The central position is marked by a cross, b. and d. Contour map of the CO(2-l) distribution and corresponding 
velocity field. The levels are the same as for the CO(1-0). 



1500 to 1960 km s" 1 at 32" SE and 72" NW, is thus 
smaller than the range of velocities in the ionized disk. 

Shaya et al. (1996) used HST V and I band images 
of the central region (a field of view of 70") to analyze 
the three-dimensional distribution of the dust. They de- 
rived a small amount of extinction toward the very center, 
Ay < 0.4 mag and Aj < 0.2 mag, and Ay ~1.5 in the 
regions of maximum extinction. The amount of molecular 
gas that we infer from the CO observations correlates with 
the amount of extinction: we detect only ^10 8 M Q in the 
center and about twice as much in the dusty NW region. 

Th e radio jet observed by Geldzahlcr & Fomalont 
( |l984|) could be deflected by interstellar gas, especially in 
the NW. Figure 11a shows a superposition of our CO(2-l) 



map and the radio jet at A20 cm. The direction of the cen- 
tral jet changes abruptly in the NW just below the peak in 
the CO distribution. There is a hint of a similar bending 
near the SE peak in the CO. 

The most luminous galaxies in the far-infrared are 
mergers. NGC 1316 was detected by IRAS in the four 
bands but its far-infrared luminosity is moderate (~2 x 10 9 
Lq). The Lfir to M{S.2) ratio, often used as an indicator 
of the efficiency at which molecular gas is converted into 
new stars, is ~5. This is lower than the average value of 
~15 that we calculated for th e sample of ellipticals ob- 
served by Wiklind et al. ( 1995 ) (after exclusion of the two 
galaxies with the highest values of Lfir/M{S-2) ratios, an 
SO and a galaxy larger than the beam size, for which the 
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H2 mass must have been underestimated). Nevertheless, 
this may not mean that the star formation efficiency of 
NGC 1316 is lower than in other ellipticals since it is 
known that some of the far-infrared emission is not re- 
lated to massive ionizing stars, but to dust heated by the 
overall stellar radiation field. Moreover, the Galactic CO- 
to-H2 conversion factor may not apply to ellipticals or it 
may vary from one elliptical to another. For example, an 
elliptical galaxy that has accreted a gas-rich but metal- 
poor companion could be CO-poor and its H2 mass would 
be underestimated. 
• The outer Hn regions 

It is interesting that SH2, despite being so much 
smaller than EELR, posses ses ab out twice as much atomic 
gas. Mackie & Fabbiano ( 1998 ) give an Ha luminosity 



are NGC 5128, NGC 3656, and Arp 230, all of which have 
disks that are viewed nearly edge-on (Schiminovich et al. 

Sofue et al. 



Ln a = 5.46xl0 5 L for the EELR, and a mass of ionized 
gas Mhii=5x 10 3 M (when corrected for our adopted dis- 
tance). If the EELR were photoionized by hot stars, this 
quantity of ionized gas could be maintained by ~ 40 O- 
type stars. However, Mackie and Fabbiano argued that 
there is no evidence of a population of young stars associ- 
ated with the EELR and suggested that the ionized gas is 
produced by shocks. Although Mackie & Fabbiano do not 
quote an ionized gas mass for SH2, it appears from their 
Ha image, which has been smoothed to a resolution of 
4'.'3, that the Ha flux is of the order of 3 x 10" 16 ergs s _1 
cm~ 2 , integrated over the l'/5 x 3'.'5 size reported by S80. 
At our adopted distance, such recombination emission can 
be maintained by a single O-type star. For comparison, 
the giant Hn region NGC 604 in the Local Group galaxy 
M 33 contains 3 x 10 6 M Q of ionized gas (Churchwell & 



Goss 1999] ) and approximately 200 O-type stars (Gonzalez 
Delgado & Perez |2000[). Although the Ha emission of SH2 



may signify recent star-forming activity, the population of 
massive stars is not large in comparison with the mass of 
neutral gas in the vicinity, Mm ~ 2 x 10 7 M Q . 

5.2. The absence of H 1 emission and absorption in the 
central spheroid 

In the central region of NGC 1316 the Hi A21 cm line is 
not clearly detected in emission nor is there any noticeable 
absorption against the continuum emission. Our limit on 
Hi emission in the spheroid of NGC 1316 could represent 
a mass of H as high as 10 8 Mg. The limit on absorption 
against the central continuum source corresponds to a col- 
umn density of H < 1 x 10 18 T S cm~ 2 , implying that only 
gas with spin temperature T s < 100 K could reasonably 
have been detected. The upper limit on Hi emission may 
at first sight seem surprising in view of the large amounts 
of dust and molecular gas in the spheroid. Fornax A is not 
unique in this respect: other well known elliptical galax- 
ies and merger remnants show CO and dust but no de- 
tectable H in their c enter s (e.g. NGC 7252, Wang et al. 
1992 ; Hibbard et al. 1994 ). In contrast there are smaller 
elliptical galaxies or merger remnants that have dust, CO 
and plenty of H all the way in to their centers. Examples 



1994; Eckart et al. 1990; Balcclls et al 



1993; Schiminovich et al. 2001; Galletta et al. 1997). The 
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H-deficient and H-rich systems also differ in the size of 
the spheroid and presence of extended X-ray emission. 
However, the detectability of Hi absorption may also de- 
pend on whether a gaseous disk is present and if so how 
it is inclined to the observer's line of sight. 

The CO emission from NGC 1316 probably arises in 
the surfaces of molecular clouds. Any neutral atomic gas 
associated with such molecular clouds may cover too little 
of the central continuum source to be detectable in ab- 
sorption in the A21 cm line at the angular resolution and 
sensitivity of the available measurements. The high pres- 
sure of the hot X-ray-emitting gas and the shocks associ- 
ated with the merger events may help confine the neutral 
gas in small clouds, while the X-radiation keeps the more 
dilute gas ionized. 

Finally, it is interesting to consider briefly how the spin 
temperature of H in the central region of NGC 1316 might 
be affected by the extreme environment there. There are 
radiative processes that can boost the spin temperature 
of low-density H: absorption and stimulated emission in- 
duced by the continuum radiation in the 21 cm line itself 
(cf. Bahcall & Ekers [1969 ), and absorption and fluores- 
cence in the Hi La line (Deguchi & Watson 1985). These 
processes will not affect the interpretation of optically thin 
Hi emission lines, but they could make absorption even less 
readily detectable. Although Fornax A is a bright contin- 
uum source at A21 cm, the integrated flux density of 0.24 
Jy over the central 10" x 30" (Ekers et al. 1983, Geldzahler 
& Fomalont 1984) corresponds to an average continuum 
brightness temperature of 565 K and yields an absorption 
rate F = 0—>loip^7x 10~ n s -1 . The brightness 
temperature in the extended central jets is lower, ~ 150 
K, which might be a lower limit on the spin temperature 
of dilute gas in the center. Pumping in the Hi La line will 
have a larger effect on the spin temperature of Hi . The 
average surface brightness of the Hi La line can be esti- 
mated from the observed surface brightness in the Ha line, 
I~4x 10~ 17 erg s" 1 cm~ 2 arcsec" 2 (Kim et al. 1998), 
with the assumption that a fraction 77 of the Balmer-a 
emissions are followed by emission in La. The correspond- 
ing pumping rate F = — > 1 is ~ 10~ 7 ?7 s" 1 , averaged 
over the central arcmin of the galaxy. Since r\ is likely to 
be of the order of unity, this excitation rate is much higher 
than the mean rate of continuum pumping in the 21 cm 
line and is high enough to dominate over spin-changing 
collisions if the density of the Hi -emitting gas is Ss 10 2 
cm -3 . When La pumping controls the spin temperature, 
this t emperature is probably ~ 10 4 K (Deguchi & Watson 
1985 ). Thus the La pumping could make the Hi line un- 
detectable in absorption in the center of NGC 1316 by en- 
suring that its spin temperature remains high even where 
the kinetic temperature might be much lower than 10 4 K. 
The true magnitude of the La pumping effect is uncertain, 
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of course, because the distributions of dilute neutral gas 
and of Lev photons are not known. 

There is a large amount of hot, ionized, X-ray-emitting 
gas in the central kpc of NGC 1316, which arises at the 
expense of neutral gas that might otherwise exist there. 
In the future, it will be valuable to obtain Hi observations 
of greater sensitivity and higher resolution to try to trace 
the neutral gas and its interactions with its environment 
in more detail. 



5.3. Timescales 

The complex structure of NGC 1316 (five tails or loops 
of varying morphology) argues for a rich dynamical his- 
tory and probably a succession of merging events. The 
presence of shells around ellipticals and minor-axis dusty 
disks can be interpreted in two different ways: as tracers of 
the accretion of a small compan ion o nto an already exist- 
ing elliptical galaxy (e.g. Quinn 1984 ), or as the returning 
of tidal material during th e me rger of t wo dis k galaxies 
(e.g., Hernquist & Spergel |l992| Balcells |l997|) . There is 
observational evidence that NGC 1316 formed in a major 
merger ^3 Gyr ago and has been accreting material more 
recently, as discussed below. 



5.3.1. A major merger ~3 Gyr ago 

NGC 1316 has a relaxed morphological appearance within 
3', which indicates that, if the galaxy is the result of a 
major merger between two disk galaxies, the event must 
have occurred more than ^2xl0 8 years ago. 

The population of globular clusters constrains the age 
of the system. Among 24 star clusters associated with 
NGC 1316, Goudfrooij et al. ( |200l| ) identified four ex- 
ceptionally luminous ones (My < 12.3) with near-solar 
metallicities, whose age of 3.0 ± 0.5 Gyr could be deter- 
mined from comparison with single-burst population mod- 
els. Those intermediate-age clusters would have formed 
during a merger of two disk galaxies that resulted in the 
formation of the elliptical whereas the older (^14 Gyr old) 
metal-poor clusters are likely to have belonged to the pro- 
genitor disk galaxies. As a remnant of the merger of two 
disk galaxies, NGC 1316 would be at a more advanced 
stage than for instance NCC 7252, which is believed to 
be 0.5-1 Gyr old and still displays extended tidal tails 



(Schweizer & Seitzer |1998[ ). 

Optical spectra of the inner 3" of NGC 1316 have re- 
vealed a stellar population of age ~2 Gyr (Kuntschner 
2000). Those stars would therefore have formed from the 
molecular gas that fell toward the center following the 
merger event and would have contributed to the deple- 
tion of molecular gas observed in the central region today. 
Stron g nuclear H/3 a bsorp tion has been found (Schweizer 
1981 ; Bosma et al. 1985|) . It might indicate an impor- 



The X-ray luminosity increases as a result of gas be- 
ing shock heated or evaporated by electrons in the already 
existing hot plasma (e.g., Goudfrooij & Trinchicri 1998[ ) 
and mass loss from evolving giant stars and supernovae la 
(e.g ., Bri ghcnti & Mathews |l999| ). Using ROSAT, Kim ct 
al. ( |1998| ) detected - 5 x 10 8 M Q (at D = 18.6 Mpc) of 
hot (~10 7 K) gas, which is small compared to the amount 
of hot gas found in some early-type galaxies of compara- 
ble optical luminosity. However, Kim et al. (1998) point 
out that the amount of hot gas is comparable with what 
would be expected from stellar evolution in ~1 Gyr. A 
comparable amount of hot gas is coincident with the tidal 
tails in the outskirts of NGC 1316 (Mackie & Fabbiano 



1998) and may be due to more recent events. 



5.3.2. More recent infalls 

The presence of loops and ripples as well as of a disk of ion- 
ized gas is strong evidence in favor of the infall of a small 
galaxy onto NGC 1316 less than 1 Gyr ago (Schweizer 



1980). In particular, the giant loop called Li by Schweizer 



is similar to features seen in simulations of shell galaxies 
(e.g., Quinn 1984, Hernquist & Quinn 1988) and could be 



interpreted as due to the infall of a dynamically colder 
galaxy ^0.5 Gyr ago. Mackie & Fabbiano (1998) have ar- 



tant component of A-type stars which are part of a post- 
starburst population. 



gued from the temperature of the hot gas detected in the 
Li region that Li could be the remnant of a small gas-rich 
galaxy on a high- velocity encounter (~380 km s _1 ) with 
NGC 1316, and its origin would therefore be different from 
that of the other loops which formed in the aftermath of 
the major, older merger. The extended emission-line re- 
gion EELR is found in the area of highest surface bright- 
ness of Li. The fact that we don't detect any CO from the 
EELR is not surprising if the infalling galaxy was metal- 
poor. The low level of detectable Hi emission from that 
region may be due to the fact that the most of the atomic 
hydrogen has been stripped, or photo-ionized by the newly 
formed OB associations. 

The NW concentration could as well be a recently 
accreted gas-rich galaxy. It has a linear dimension of 
about 2.5 kpc and contains about 2.2xl0 8 Mq of molec- 
ular hydrogen. This H2 mass is comparable to that of 
NGC 1317, the small (~15 kpc diameter) spiral compan- 
ion of NGC 1316 (Horellou et al. |1995[ ). This gas concen- 
tration is located 40" (3.5 kpc) from the nucleus, slightly 
to the north of the minor axis. The velocity of the molec- 
ular gas, v ^1850 km s _1 , is ~70 km s _1 higher than 
the systemic velocity given by Arnaboldi et al. (1998) and 
agrees exactly with the stellar velocity derived from ab- 
sorption line spectra in that region (figure 3 of Arnaboldi 
et al. 1998). Interestingly, beyond 30" from the nucleus 
along the minor axis, the stellar radial velocities increase 
on both sides, producing a U-shaped velocity curve, and 
the velocity dispersion profile is irregular. Arnaboldi et al. 
argue that this kinematical feature is caused by infalling 
material within one dynamical time-scale at 40" (i.e. less 
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a. CO(2-l)map 
+ 20 cm radio jet 



b. CO(2-l)map 

+ sketch of dust lanes 



c. HST B-l image 



Fig. 11. a. CO(2-l) contour map overlaid on the 20 cm radio jet mapped by Geldzahler & Fomalont (1984). b. 
CO(2-l) contour map overlaid on Schweizer (1980)'s sketch of the dust l anes (dots), c. HSTB-V image of the central 
region of NGC 1316 (press release STScI-PR99-06, Grillmair et al. |l999|) . 



than 3 x 10 7 years ago) 
pected to last longer. 



since such features are not ex- 



5.4. Gaseous shells? 

It is interesting to compare NGC 1316 with NGC 5128 
(Centaurus A). NGC 5128 lacks tidal tails but it has a 
more prominent dust lane and a better defined system of 
shells in both of which atomi c and molecular gas have 
been detected (Eckart et al. 1990|; Schiminovich et al. 



199J; Charmandaris et al. |2000|) . NGC 5128 has com- 
parable amounts of Hi and H2 in the region of the dust 
lane (~ 3.5 x 10 8 M ) whereas NGC 1316 has at least 
four times more molecular gas than Hi . This difference in 
Hi content may be due to the fact that the Hi in NGC 1316 
is more affected by the X-ray emission, since NGC 1316 is 
brighter in X-ray than NGC 5128 (factor of 8 in the 0.1-2.4 
keV band, ROSAT observations, and factor of 1.4 in the 
0.2-4 keV band, Einstein IPC observations). Both galax- 
ies have a double-lobe radio source. The radio distribution, 
however, has a different orientation: in Centaurus A it is 
perpendicular to the dust lane whereas in Fornax A it is 
not. Nuclear effects may therefore affect the gas and dust 
properties more clearly in Fornax A than in Centaurus A. 
In the outer parts of NGC 5128, Hi is found in the shells, 
with CO detected in the highest column density regions 
with M(Hi )/M(H 2 )~l. In NGC 1316, we could not detect 
any CO from the bright outer H11 regions, but our limit on 
the H2 mass is comparable to the Hi mass. The possible 
chain of Hi clumps in the outer part of NGC 1316 could 
be the highest peaks in a continuous shell. More sensitive 
observations are required to confirm this. 



The HST image of NGC 1316 displays beautifully the 
distribution of the dust (Fig. lie). Northwest of the nu- 
cleus it roughly has a T-shape with a radial extent of about 
40" and an azimuthal extent of almost 90°. CO is detected 
in the area of both the radial and the azimuthal dust 



patches. Grillmair et al. (1999) point out the similarity 



of the radial structure to the photo-evaporating columns 



in the Eagle Nebula M16 (Hester et al. 1996), though on 



a much larger scale, and suggest that it may be due to 
Rayleigh- Taylor instabilities or subject to erosion by par- 
ticles or radiation from the nucleus of the galaxy. The 
azimuthal feature may be a dusty shell. The azimuthal 
configuration is expected from recently accreted mate- 
rial oscillating between two apocenters and forming shells 
(e.g. Hernquist & Quinn 1988). Gaseous shells could form 
in the same phase-wrapping process if the gas is clumpy 
and not very dissipative (Charmandaris & Combes 2000 ). 
A detailed dynamical model of the gaseous features in 
NGC 1316 based on the accretion of a small companion 
will be the subject of another paper. 

6. Summary and conclusions 

Our new observations of molecular gas (as traced by 
12 CO(1-0) and 12 CO(2-l) line emission) toward the pe- 
culiar elliptical galaxy NGC 1316 and a new analysis of 
Hi data have shown that: 



1. 



About 5 x10 s Mq of molecular gas is located in the 
inner 2' of the galaxy and is mainly associated with 
the dust patches along the minor axis (~1.8xl0 8 M Q 
in the southeast, ~10 8 M Q in the center and ~2.2xl0 8 
M Q in the northwest). Those mass estimates are based 
on a Galactic CO/H 2 conversion factor. 
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2. No neutral atomic gas found was found in the spheroid 
of NGC 1316, down to a limit of -10 8 M©. Hi was de- 
tected at four locations in the outer part of NGC 1316. 

3. Two of those locations coincide with bright Hn regions 
in a tidal loop. About 2xl0 7 M Q of Hi was found 
toward the giant Hn region identified by Schweizer 
( 1980Q and located 6'7 (or 36.2 kpc) from the nucleus 
of NGC 1316. No CO could be detected. About 10 7 
M Q of Hi was detected toward the extended emission 
line region studied by Mackie & Fabbiano ( 1998J ) . 

4. The NW molecular gas concentration could be due to 
recent infall of material onto NGC 1316. It could be 
a small gas-rich galaxy or the trace of a gaseous shell 
formed though radial oscillations of infalling material. 

5. The two spiral companions of NGC 1316, NGC 1317 
to the north and NGC 1310 to the east, are lo- 
cated within the primary beam of the VLA and their 
Hi distribution and kinematics could be mapped. It 
was found that NGC 1317 has an unusually small 
Hi disk that may have been affected by ram-pressure 
stripping if it passed through the X-ray emitting re- 
gion of NGC 1316. The Hi distribution of NGC 1310 
is asymmetric. 

6. There is a hint that the central radio jet of Fornax A is 
bent owing to interaction with the large concentrations 
of interstellar gas seen in CO. 
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